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Executive  Summary 

Over  the  approximately  2  Vz  year  period  of  the  above  AFOSR  grant,  substantial  progress  has  been 
made  in  the  further  development  of  Ballistic  Electron  Emission  Microscopy  (BEEM)  as  a 
quantitative  microscopic  and  spectroscopic  tool  for  the  characterization  of  buried  quantum  objects 
and  electronic  defects  in  semiconductors.  During  the  period  of  this  grant,  we  have  done  the 
following: 

i.  Developed  a  theoretical  model  for  the  simulation  of  BEEM  spectra.  The  model  is  compared 
with  the  experimentally  obtained  Second  Derivative  (SD)  spectra  for  the  prototypical  case  of 
AlxGai_xAs/GaAs  heterostructures.  Quantitative  numerical  agreement  is  obtained. 

ii.  We  have  extended  the  modeling  with  Monte  Carlo  simulations  to  get  detailed  understanding 
of  BEEM  spectra  as  a  function  of  depth  below  the  surface. 

iii.  BEEM  has  been  used  for  several  spectroscopic  studies  on  novel  materials  systems.  These 
include  observation  of  the  splitting  of  the  L  band  in  GalnP  due  to  ordering;  the  “local” 
conduction  band  offset  of  GaSb  quantum  ‘dots’  grown  on  GaAs;  study  of  local  electrical 
transport  within  the  network  of  dislocations  in  GaN;  and  finally,  the  measurement  of 
heterojunction  band  offsets  in  several  planar  heterostructures  based  on  arsenides,  phosphides 
and  antimonides. 

The  work  reported  here  was  performed  primarily  by  M.  Rubin  (PhD  ’98),  S.  Bhargava  (PhD  ’98), 
E.Y.  Lee  (Assistant  Research  Engineerg),  M.  Kozhevnikov  (Research  Associate),  C.  Zhang 
(Graduate  Student),  and  H.H.Lee  (Graduate  Student). 
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I.  Introduction 


Ballistic  Electron  Emission  Microscopy  (BEEM)  is  a  three-terminal  modification  of  Scanning 
Tunneling  Microscopy  (STM)  which  has  been  widely  used  for  the  study  of  Metal-semiconductor 
(m-s)  interfaces.1,2,3  More  recently 4  the  Principal  Investigator  of  this  grant  and  his  group  have 
extensively  used  this  technique  to  study  semiconductor  heterojunctions  buried  substantially  below 
the  m-s  interface. 

Over  the  last  three  years  we  have  made  substantial  progress  in  developing  quantitative  theoretical 
models  for  BEEM  and  compared  them  with  our  experiments,  as  well  as  gained  understanding  of  the 
scattering  processes  inside  the  semiconductor.  Finally  several  new  materials  systems  have  been 
studied.  This  final  Technical  Report  provides  a  summary  of  this  work. 


II.  Scattering  Models  and  Experiment 

The  original  theoretical  description  of  BEEM  was  based  on  a  planar  tunneling  model  and  assumed 
transverse  momentum  conservation  2  at  the  m-s  interface.  However,  until  recently  it  was  not 
possible  to  quantitatively  describe  the  BEEM  spectra  anisotropy  or  distribution  between  different 
conduction  channels  for  non-epitaxial  systems  such  as  Au/GaAs  and  Au/Si. 

To  analyze  quantitatively  the  experimental  data,  the  m-s  interface-induced  scattering  (MSIS)  model 
was  proposed  by  us.5  The  authors  showed  that  the  experimental  data  for  Au/Si  and  Au/GaAs 
systems  could  be  fitted  only  by  taking  into  account  of  the  strong  carrier  scattering  at  the  m-s 
interface.  Recently,  to  quantitatively  analyze  the  electron  scattering  at  the  m-s  interface,  this  model 
was  modified  to  include  anisotropy  of  the  electron  effective  mass,  the  energy  dependence  of  the 
electron  mfp  in  the  metal  base,  and  finite  temperature.6  In  this  model,  the  conduction  process  is 
described  by  carrier  multivalley  transport  with  corrections  due  to  the  scattering  at  the  m-s  interface, 
which  depends  on  the  m-s  scattering  probability.  When  interface  scattering  is  considered,  the 
BEEM  current  is  given  by  6 


lc  =  eA Y  F[ (ft )TV‘  (ft ,  ft )D(ft ) 


X 


-  ^  pp  f>f  y '  (ft ,  ftb  ^  (pb  )+jJpfi£fT '  & .  ftb  )t  h(ftb) 


(2) 


where  P^b  is  the  probability  for  the  electron  with  the  wavevector  kb  to  be  scattered  into  the  state 

labeled  by  the  wavevector  ft .  The  first  and  second  terms  in  brackets  correspond  to  electrons  that 
are  not  scattered  and  that  are  scattered,  respectively. 


4 


Au/Si 


Fig.  1.  Comparison  of  the  calculation  in  the  strong  scattering  limit  for  Au/Si(100)  (dashed  line)  and  Au/Si(lll)  (solid 
line)  with  the  experimental  results  for  Au/Si.  These  measurements  were  made  on  Au/Si(l  1 1)  but  measurements  for 
Au/Si(100)  are  very  similar  to  those  for  Au/Si(lll).  The  inset  compares  the  calculation  in  the  weak  scattering  limit  to 
the  same  experimental  results  (from  D.  L.  Smith,  E.  Y.  Lee,  and  V.  Narayanamurti,  Phys.  Rev.  Lett.  80,  2433  (1998)). 


In  this  model,  if  scattering  at  the  m-s  interface  is  strong,  the  injected  electron  flux  distribution  is 
redistributed  and  valleys  with  zero  interface  transverse  wavevector  at  their  energy  minimum  are  not 
preferentially  weighted.  Instead  the  weighting  goes  like  the  density  of  final  states  for  the  scattering 
process.  Fig.  1  compares  the  calculated  BEEM  current  in  the  strong  and  in  the  weak  scattering 
limits  for  Au/Si  with  the  experimental  results.  One  can  see  that  the  experimental  data  can  be 
described  adequately  without  fitting  parameters  only  in  the  strong  scattering  limit.  In  BEEM  on 
Au/GaAs,  the  calculated  T  channel  contribution  is  much  larger  than  the  L  channel  contribution  to 
the  collector  current  in  the  weak  scattering  limit,  whereas  the  L  channel  contribution  is  much  larger 
than  the  F  channel  contribution  in  the  strong  scattering  limit  (see  Fig.  2).  Scaled  threshold  shapes 
are  similar  in  the  weak  and  strong  scattering  limits  for  conduction  band  valleys  that  make  a  circular 
projection  on  the  interface  plane,  such  as  F  and  X  valleys.  For  conduction  band  valleys  that  make  a 
strongly  anisotropic  projection  on  the  interface  plane,  such  as  the  L  valleys  of  GaAs(lOO),  the 
calculated  L  channel  currents  have  significantly  different  threshold  shapes  in  the  two  limits,  as 
shown  in  Fig.  3.  When  the  calculations  are  done  in  a  spherical  approximation,  the  threshold  shapes 
are  similar  in  the  two  limits  and  look  like  the  result  of  the  more  complete  anisotropic  model  in  the 
strong  scattering  limit. 
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Fig.  2.  Calculated  ratio  of  collector  to  tunneling  current  as  a  function  of  bias  for  the  three  conduction  band  channels  of 
Au/GaAs(100)  in  the  weak  scattering  limit  (upper  panel)  and  in  the  strong  scattering  limit  (lower  panel)  (from  D.  L. 
Smith,  M.  Kozhevnikov,  E.Y.  Lee  and  V.  Narayanamurti,  Phys.  Rev.  B61,  13914  (2000)). 


Bias  (V) 


Fig.  3.  Scaled  ratio  of  collector  to  tunneling  current  (upper  panel)  and  scaled  second  voltage  derivative  of  the  ratio  of 
collector  to  tunneling  current  (lower  panel)  as  a  function  of  bias  for  the  three  conduction  band  channels  of 
Au/GaAs(100).  The  scaling  factors  for  all  the  valleys  were  chosen  so  that  the  currents  for  two  limits  overlap  well  in  the 
near  threshold  region.  The  strong  scattering  limit  results  (solid  lines)  are  as  calculated  and  the  weak  scattering  limit 
results  (dashed  lines)  are  scaled  by:  Z  channel,  divide  by  31;  L  channel,  multiply  by  4.5;  X  channel,  divide  by  1.7  (from 
D.  L.  Smith,  M.  Kozhevnikov,  E.Y.  Lee  and  V.  Narayanamurti,  Phys.  Rev.  B61,  13914  (2000)). 


6 


Recently,  the  electron  scattering  at  the  metal-semiconductor  interface  was  studied  experimentally 
by  us  in  the  Au/GaAs/AlxGai.xAs  system  by  using  second  derivative  (SD)  BEEM  spectroscopy.7 
The  SD-BEEM  spectrum  is  proportional  to  the  heterostructure  transmission  coefficient,  and  is 
particularly  useful  both  for  separating  the  contributions  from  various  conduction  band  channels  and 
for  analyzing  the  threshold  shapes.  The  analysis  of  the  spectral  shape  is  more  effective  than  just  the 
BEEM  current  magnitude  because  the  absolute  value  of  the  BEEM  current  could  be  scaled  by  the 
fluctuations  in  both  the  metal  layer  and  semiconductor  epilayer  thickness.  In  GaAs,  conduction 
bands  with  both  zero  and  large  interface  parallel  wavevector  coexist  in  the  same  material.  Thus 
direct  and  indirect  transport  channels  can  be  investigated  in  the  same  sample.  The  observed  SD- 
BEEM  spectral  shape  shows  obvious  deviation  from  the  calculations  that  assume  transverse 
momentum  conservation.  In  Fig.  4,  we  compare  the  calculated  (dashed  lines)  collector  current  and 
its  second  voltage  derivatives  with  the  corresponding  measured  quantities  (solid  lines)  for  a 
Au/GaAs  sample  at  room  temperature.  The  SD  spectrum  was  extracted  from  the  experimental 
BEEM  spectrum  in  Fig.  4(a)  by  numerical  differentiation  with  a  10  meV  window,  as  shown  in  Fig. 
4(b). 
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Fig.4  The  room-temperature  BEEM  (a)  and  SD-BEEM  (b)  spectra  of  the  1  pm  undoped  GaAs  layer  grown  on  n+  GaAs 
substrate.  The  MSIS  model  calculations  (dashed  lines)  are  also  presented  for  three  values  of  SP,  the  electron  scattering 
probability  at  the  m-s  interface.  The  inset  in  (a)  shows  a  schematic  band  diagram  of  the  samples  under  study  (from  M. 
Kozhevnikov,  V.  Narayanamurti,  C.  Zheng,  Y.-J.  Chiu,  and  D.  L.  Smith,  Phys.  Rev.  Lett.  82,  3677  (1999)). 


The  probability  of  electron  scattering  at  the  m-s  interface  (SP)  was  adjusted  to  fit  the  SD-BEEM 
spectra.  The  SD-BEEM  spectra  clearly  separate  the  contributions  from  T  and  L  electrons;  thus,  the 
weighting  of  the  relative  T  and  L  channel  contribution  is  a  sensitive  test  of  the  model.  The  best  fit 
to  the  SD-BEEM  spectrum  gives  an  85%  probability  for  electron  scattering  at  the  m-s  interface.  For 
comparison,  theoretical  curves  for  the  cases  without  scattering  and  with  100%  scattering  are  also 
shown.  The  measured  relative  contribution  of  the  L-electrons  is  different  than  that  calculated  for 
the  L  channel  in  the  weak  scattering  limit.  Specifically,  the  L-electron  collector  current  is  found  to 
be  the  strongest  one,  whereas,  assuming  transverse  momentum  conservation,  the  contribution  of  the 
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L  band  (which  does  not  project  to  the  zone  center  of  the  interface  Brillouin  zone)  can  give  only  a 
small  BEEM  current  onset.  The  T  and  L  channel  contributions  are  included  in  the  calculation  but 
the  X  channel  contribution  is  not  included.  In  this  experimental  structure,  the  absence  of  the  X 
conduction  minimum  contribution  at  room  temperature  is  due  to  the  image  potential-induced  strong 
electron  scattering  in  the  spacer  between  the  metallurgical  m-s  interface  and  the  maximum  of  the 
barrier  height.  Fig.  5  presents  the  MSIS  model  fits  to  the  BEEM  and  SD-BEEM  spectra  for  a 
number  of  Au/GaAs/AlxGai.xAs  single  barrier  structures.  In  this  case,  the  best  fits  are  obtained  with 
the  scattering  probability  at  the  m-s  interface  varying  between  85%  and  92%  among  these  samples. 
This  small  variation  in  the  scattering  parameter  indicates  that  our  diode  fabrication  procedure  is 
reproducible  and  results  in  approximately  the  same  quality  of  the  m-s  interface. 


SB  structure  (G  aAs/AI^G  a,  xAs/G  aAs) 


Fig.  5.  Room-temperature  SD-BEEM  spectra  for  five  different  A1  compositions  (solid  lines).  For  clarity,  the  SD- 
BEEM  spectra  are  shifted  along  the  vertical  axis.  Thin  solid  lines  are  eye-guides  for  the  peaks  position  development. 
The  MSIS  model  calculations  are  also  presented.  The  model  calculations  show  the  separate  T  and  L  valley  contributions 
(dotted  lines)  and  their  sum  (dashed  lines)  (from  M.  Kozhevnikov,  V.  Narayanamurti,  C.  Zheng,  Y.-J.  Chiu,  and  D.  L. 
Smith,  Phys.  Rev.  Lett.  82,  3677  (1999)). 


In  summary,  through  the  combination  of  SD-BEEM  spectroscopy  and  theoretical  modeling,  which 
takes  into  account  the  parallel  momentum  non-conversation,  we  are  able  to  quantitatively  describe 
BEEM  data.  This  ability  clearly  makes  BEEM  a  more  useful  spectroscopic  and  microscopic  tool. 
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III.  BEEM  Depth  resolution 


The  hot  electron  transport  in  BEEM  experiments  can  be  considered  to  be  ballistic  only  with  certain 
restrictions.  The  real  picture  is  certainly  different  because  of  consecutive  scattering  processes  in  the 
metal  base  film,  at  the  metal/semiconductor  interface  and  in  the  semiconductor  itself.  Even  for  the 
case  of  a  perfectly  clean  device  and  extremely  low  temperatures,  the  carrier  mean  free  path  is 
usually  comparable  with  the  device  thickness.  Therefore,  the  ballistic  description  of  the  electron 
transport  is  restricted  to  layers  close  to  the  surface.  As  a  consequence,  BEEM  can  be  characterized 
by  its  depth  resolution. 
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Fig.  6.  For  the  populations  of  the  electrons  in  GaAs  resulting  from  a  BEEM  experiment  of  Au-GaAs,  the  energy  and 
angular  distributions  are  shown.  The  y  axis  is  the  cosine,  denoted  by  p,  of  the  angles  of  the  velocities  of  the  electrons 
with  respect  to  the  interface  normal  and  the  X  axis  is  the  kinetic  energy  E  of  the  electrons  with  respect  to  the  T- valley 
minimum.  The  condition  is  the  strong  scattering  limit  of  SLN  theory  at  300  K  and  at  a  bias  of  1.53  V.  Initially,  1000 
electrons  are  injected  at  the  Au-GaAs  interface,  (a)  The  initial  electron  distributions  of  T-  (left  panel  and  crosses),  L- 
(middle  panel  and  circles),  and  X-  (right  panel  and  diamonds)  valley  electrons  at  the  Au-GaAs  interface,  (b)  The 
distributions  at  600  A  below  the  Au-GaAs  interface  (from  E.  Y.  Lee,  Phys.  Rev.  59,  15332  (1999)). 


Monte  Carlo  simulations  8  of  electron  transport  in  the  T,  L  and  X  valleys  of  GaAs  were  performed 
using  the  MSIS  model.5  With  increasing  object  depth  beneath  the  Au/GaAs  interface,  the  Monte 
Carlo  simulations  predict  (1)  significant  cooling  of  hot  electrons,  on  the  order  of  ~3meV/nm  at  1.5 
eV,  and  (2)  significant  redistribution  of  electrons  among  the  conduction  bands  (see  Fig.  6).  In  order 
to  show  experimentally  that  scattering  inside  the  semiconductor  structure  affects  the  BEEM  spectra, 
in  Ref.  7  we  compared  the  BEEM  spectra  for  several  pairs  of  Au/GaAs/AlGaAs  single  barrier  (SB) 
samples  with  the  same  SB  composition  but  with  different  cap  thickness  of  50  A  and  300  A.  The 
obtained  SD-BEEM  spectra  of  the  GaAs/Alo.iGao.gAs/GaAs  SB  samples  are  shown  in  Fig.  7  for 
T=85K  and  300K  (these  data  are  representative  for  all  samples  pairs). 
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Fig.  7.  SD-BEEM  spectra  for  GaAs/Alo.2Gao.sAs/GaAs  SB  samples  with  50  A  and  300  A  GaAs  cap  layers,  taken  at 
7=300  K  (solid  curves)  and  T- 85  K  (dotted  curves).  For  clarity,  the  spectra  are  shifted  along  the  vertical  axis  (adopted 
from  M.  Kozhevnikov,  V.  Narayanamurti,  C.  Zheng,  Y.-J.  Chiu,  and  D.  L.  Smith,  Phys.  Rev.  Lett.  82,  3677  (1999)). 


Since  the  expected  mfp  lengths  for  the  T,  L  and  X  electrons  at  T=300K  are  -1000  A,  -100  A  and 
-10  A,  respectively, 9  manipulations  with  the  heterostructure  thickness  in  the  range  of  100-300  A 
should  affect  mainly  the  L-electrons  contribution  to  the  BEEM  current  (the  effective  heterostructure 
thickness,  that  will  affect  the  BEEM  current  collection,  is  the  combined  thickness  of  the  GaAs  cap 
layer  and  the  SB  layer.)  Indeed,  as  one  can  see  from  Fig.  7,  the  L-electrons  contribution  for  the 
sample  with  300  A-cap  layer  (with  the  SB  thickness  of  50  A,  the  total  effective  heterostructure 
thickness  is  -350A)  is  reduced  by  factor  of  3  as  compared  with  the  50  A-cap  layer  sample  (the  total 
effective  thickness  of  -100  A),  whereas  the  T  contribution  has  remained  intact.  As  the  temperature 
decreases  from  300  K  to  85  K,  in  addition  to  the  spectrum  shift  expected  from  the  temperature 
dependence  of  the  energy  gap,  a  strong  increase  of  the  signal  is  observed  for  the  L-electrons  in  the 
SB  sample  with  the  300  A-cap  layer,  whereas  the  SD-BEEM  spectral  shape  remains  essentially  the 
same  for  the  SB  sample  with  the  50  A-cap  layer.  These  experimental  results  are  in  agreement  with 
the  expected  increase  in  the  mfp  of  the  electrons  with  the  decreasing  temperature.  As  the 
temperature  decreases  from  300  K  to  85  K,  the  calculated  mfp  near  the  energy  threshold  increases 
from  -1000  A  to  -1500  A  for  T  electrons  and  from  -100  A  to  -300  A  for  the  L  electrons.9  The  X- 
electrons  contribution  is  completely  obscured  in  the  room-temperature  BEEM  spectra  presented  in 
Fig.  7  due  to  their  cooling-induced  effective  scattering  in  the  GaAs  cap  layer.  Indeed,  the  X- 
electrons  mfp  (-10  A  at  T=300  K  and  -30  A  at  T=85  K)  is  shorter  than  the  effective  heterostructure 
thickness  for  all  samples  (>100  A).  Thus,  while  the  initial  electron  distribution  among  the 
conduction  bands  of  the  semiconductor  is  specified  by  transport  through  the  metal  base  layer  and  by 
the  m-s  interface  scattering,  further  electron  transport  is  governed  by  the  difference  in  the  electron 
mean-free-path  (mfp)  length  for  the  T,  L  and  X  electrons.  In  imaging  the  buried  semiconductor 
structures,  the  mfp  length  (in  the  case  of  GaAs-based  structures,  the  mfp  length  of  T-electrons  as  the 
longest  one)  will  determine  the  depth  resolution  of  BEEM.  It  implies  that  the  BEEM  contrast  will 
be  not  enough  to  resolve  a  semiconductor  structure  buried  deeper  than  -1000A  (mfp  of  the  T- 
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electrons).  Therefore,  as  the  heterostructure,  characterized  by  several  transport  channels,  is  buried 
deeper,  the  information  about  these  conduction  channels  is  gradually  gone,  starting  from  the 
transport  channel  with  the  shortest  electron  mfp  length. 

Lee  et  al 10  have  performed  Monte  Carlo  simulations  of  the  transport  of  electrons  injected  into  the  T 
valley  of  GaAs  for  BEEM  imaging  and  spatially  resolved  spectroscopy  of  model  quantum  dots 
(QD)  and  quantum  wires  (QW)  buried  beneath  the  Au/GaAs  interface.  For  Monte  Carlo 
simulations,  a  modified  BK  model  was  used,  with  the  buried  object  forming  a  local  barrier  for  the 
electron  transmission  (such  as  AlSb  and  GaSb  dot  grown  on  GaAs).  The  collection  plane  was 
located  at  2000A  below  the  Au/GaAs  interface,  and  the  QD  (QW)  was  modeled  as  a  thin  disk  300A 
in  diameter  (a  thin  300Awide  stripe),  which  specularly  scatters  electrons  of  all  energies.  Fig.  8 
shows  the  calculated  BEEM  current  as  a  function  of  x,  a  lateral  displacement  between  the  STM  tip 
and  the  center  of  the  buried  dot  (wire),  for  three  depths  of  the  QD  (QW)  from  the  Au/GaAs 
interface.  One  can  see  that  the  BEEM  contrast  decreases  with  the  depth  of  the  object.  At  a  depth  of 
600A,  a  QD  is  hardly  detectable  by  BEEM  whereas  the  QW  still  shows  significant  contrast,  due  to 
the  added  dimension.  In  addition,  it  was  shown  that  due  to  an  interplay  between  a  geometric 
filtering  effect  and  the  role  of  the  electron-phonon  scattering  (assuming  electron  transverse 
momentum  conservation  at  the  m-s  interface),  there  is  a  range  of  optimum  depths  for  the  sharpest 
crystal  momentum  and  energy  distribution  of  the  electrons  incident  upon  the  buried  structure. 
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Fig.  8.  The  BEEM  current  as  a  function  of  lateral  displacement  at  a  tip-to-sample  bias  of  1  V  (a)  and  1.25  V  (b)  for  a 
perfectly  reflecting  quantum  dot  300  A  in  diameter,  and  at  a 

tip-to-sample  bias  of  1  V  (c)  and  1.25  V  (d)  for  a  perfectly  reflecting  300A  wide  quantum  wire.  The  lines  are  drawn  as  a 
guide  for  the  eye.  In  all  cases,  the  solid  line  (circle)  is  for  d=100  A,  the  dashed  line  (square)  is  for  d= 300  A,  and  the 
dotted  line  (diamond)  is  for  d-600  A  (adopted  from  E.  Y.  Lee,  V.  Narayanamurti,  and  D.  L.  Smith,  Phys.  Rev.  B55, 
R 16033  (1997)). 
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As  an  example  of  BEEM  imaging  of  deeply  buried  localized  objects,  a  study  of  misfit  dislocations 
at  the  InxGai.xAs/GaAs  interface  buried  400-700  A  below  the  m-s  interface  was  conducted  in  our 
laboratory.9’1112  The  cross-hatch  misfit  dislocation  network  in  InGaAs/GaAs  was  ideally  suited  for 
BEEM  as  the  dislocation  core  beneath  the  surface  gives  rise  to  a  perturbation  at  the  surface  easily 
visible  by  AFM  or  STM.  Samples  consist  of  InGaAs  of  different  thickness  grown  on  (001)  GaAs 
substrates  and  were  designed  such  that  the  dislocations  will  exist  (In  concentrations  of  x=0. 15-0.3 
and  InGaAs  thickness  of  d=400-750A  were  examined).12  Simultaneous  STM  and  BEEM  imaging, 
shown  in  Fig.  9,  revealed  a  decrease  in  the  BEEM  current  directly  associated  with  the  cross-hatch 
pattern  in  the  STM. 


STM 


BEEM 


16A  full  scale  5  pA  full  scale 

Fig.  9.  lpxlp  STM  (left)  and  BEEM  (right)  images  of  InGaAs/GaAs  heterostructure  (x=0.2,  d=600A).  The  STM 
image  shows  the  growth  perturbation  caused  by  the  dislocation  core.  The  BEEM  image  shows  the  decrease  in  current  at 
the  presence  of  the  cross-hatch  (from  S.  Bhargava,  Ph.D.  dissertation,  University  of  California  at  Santa  Barbara,  August 
1998). 


The  minimum  in  the  BEEM  current,  through  analysis  of  line  scans,  was  found  to  be  shifted  -400  A 
from  the  maximum  of  the  STM  image,  indicating  that  the  loss  in  BEEM  current  is  not  simply  a 
surface  affect,  but  one  that  arises  from  the  dislocation  core  located  ~600A  beneath  the  surface  (see 
Fig.  10).  This  is  in  agreement  with  the  theory  that  the  dislocation  core  comes  to  the  surface  not 
along  the  growth  direction,  but  along  the  glide  plane  (111),  which  is  54.7°  off  of  the  (100)  direction. 
Monte  Carlo  calculations  were  used  to  understand  the  mechanism  responsible  for  this  decrease  12. 
The  simulation  of  the  charged  model  fit  closest  to  the  experimental  data  suggesting  that 
backscattering  off  of  line  charges  to  the  metal  base  is  responsible  for  the  loss  in  BEEM  current. 


12 


Distance  (A) 


Fig.  10.  STM  (left)  and  BEEM  (right)  linescanes  taken  over  dislocation  images.  The  STM  linescanes  show  a  Gaussian 
peak  while  the  BEEM  linescanes  show  a  drop  in  the  BEEM  current  shifted  from  the  peak  in  the  STM  image  (from  S. 
Bhargava,  Ph.D.  dissertation,  University  of  California  at  Santa  Barbara,  August  1998). 


IV.  GalnP 

In  the  past  decade,  spontaneous  CuPt  ordering  of  many  m-V  alloys  has  been  widely  observed  in 
vapor  phase  growth  on  (001)  substrates.  In  GaInP2,  extensive  theoretical 13,14  and  experimental 15 
work  has  been  carried  out  to  study  the  effect  of  the  ordering-induced  reduction  of  the  crystal 
symmetry  on  the  structural,  optical  and  transport  properties  of  the  ordered  material.  Ordering 
induced  changes  in  the  band  structure  of  GaInP2  are  important  for  its  application  in  advanced  solar 
cells  with  very  high  conversion  efficiency,  for  potential  band-gap  engineering  and  for  fundamental 
studies  of  atomic  ordering. 

Perfectly  ordered  material  has  not  been  observed.  The  electronic  states  of  a  partially  ordered 
structure  can  be  interpolated  from  those  of  the  totally  disordered  and  the  perfectly  ordered 
semiconductor.14  The  degree  of  ordering  is  usually  described  by  the  ordering  parameter  r),  where 
Eg(r|)=Eg(r|=0)-  r|2AEg(r|=l),  with  r)=0  and  r|=l  to  describe  perfectly  disordered  and  perfectly 
ordered  material,  respectively.  Eg(r|=0),  the  bandgap  of  disordered  GalnP,  is  2.01  eV  at  low 
temperature  16  and  AEg(ri=l),  the  maximum  bandgap  reduction  for  perfectly  ordered  GalnP,  is  0.47 
eV.17  In  the  most  ordered  GaInP2  structures,  the  highest  ordering  parameter  7]  -  0.6, 18  and  the 
degree  of  ordering  is  nonuniform  on  the  local  scale. 

To  characterize  the  structural  and  electronic  properties  of  GaInP2,  we  have  applied  the  BEEM 
technique.19, 20  The  undoped  GaAs/GaInP2  structures  on  n+-  and  p+-  GaAs  substrates  were  grown 
by  metal  organic  chemical  vapor  deposition  (MOCVD),  at  Tg  =  650°C.  The  structures  consist  of  a 
500  A  undoped  GaAs  buffer  layer,  a  1  pm  GaInP2  and  a  50A  GaAs  cap  layer.  A  detailed  analysis  is 
presented  here  for  the  GaInP2  samples  grown  on  a  (511)  GaAs  substrate  and  on  a  (001)  GaAs 
substrate  misoriented  6°  toward  [1 1 1]B-  From  the  low-temperature  photoluminescence,  the  bandgap 
energy  (after  correction  on  the  exciton  binding  energy)  is  -2.00  eV  for  GaInP2  grown  on  a  (51 1) 
GaAs  substrate  and  -1.88  eV  for  GaInP2  grown  on  a  6°[1 1  1]b-(001)  GaAs  substrate.  Therefore, 
we  conclude  that  the  GaInP2  layer  grown  on  a  (511)  GaAs  substrate  is  highly  disordered  (hereafter, 
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disordered  GaMy  and  the  GaInP2  grown  on  a  (001)  GaAs  substrate  misoriented  6°  toward  [111]b 
is  highly  ordered,  rj  ~  0.5  (hereafter,  ordered  GalnPa).14 

Fig.  10  shows  room-temperature  1  pm  x  1  pm  STM  images  of  the  disordered  (a)  and  ordered  (b) 
GaInP2  layer  grown  on  n+-GaAs  substrates.  One  can  see  that  the  surface  of  ordered  material  forms 
[110]  -oriented  steps.  For  ordered  GaInP2grown  on  misoriented  substrates,  the  [110]  steps  are 
usually  observed  to  form  only  single  CuPt  variant.21,22  For  disordered  GaInP2,  the  surface  is  found 
to  be  much  flatter  than  that  of  ordered  GaInP2.  The  complementary  BEEM  images  are  also  shown 
in  Fig.  10.  The  observed  contrast  of  the  BEEM  image  for  the  ordered  GaInP2  sample  is  in  direct 
correlation  with  the  surface  morphology,  indicating  high  sensitivity  and  high  spatial  resolution  of 
the  BEEM  technique. 


(a)  (b) 


Fig.  10.  Room-temperature  1  pm  x  1  pm  STM  image  (top)  and  BEEM  image  (bottom)  of  GaAs/GaInP2  layer  grown  by 
MOCVD  on  (a)  (511)  n-GaAs  substrate,  (b)  on  n-GaAs  substrate  misoriented  by  6°  to  [11 1]B-  A  tip  bias  is  -1.7  V  and  a 
tunnel  current  is  4  nA  (from  M.  Kozhevnikov,  V.  Narayanamurti,  A.  Mascarenhas,  Y.  Zhang,  J.M.  Olson,  and  D.  L. 
Smith,  Appl.  Phys.  Lett  75, 1 128  (1999)). 


To  study  the  heterostructure  transmission  coefficient,  we  analyze  the  second  voltage  derivative  (SD) 
of  the  BEEM  spectra  rather  than  the  original  BEEM  spectra.20  Fig.  1 1  shows  the  SD-BEEM  spectra 
obtained  from  the  experimental  BEEM  data  by  numerical  differentiation  with  a  10  meV  window. 
We  associate  two  clearly  pronounced  features  in  the  77K  SD-BEEM  spectra  of  disordered  GaInP2 
with  the  /’and  L  conduction  minima  contribution  in  GaInP2.  Theoretical  fits  to  the  SD-BEEM 
spectrum  for  the  disordered  GaInP2  sample,  using  the  MSIS  model,5’7  are  shown  in  Fig.  1 1  by  the 
dashed  line.  The  MSIS  model-fit  describes  the  experimental  BEEM  spectrum  reasonably  well, 
giving  -90%  probability  of  the  electron  scattering  at  the  m-s  interface,  similar  to  our  previous 
results  for  GaAs/AlGaAs  structures.7  The  absence  of  a  contribution  from  the  X  conduction 
minimum  is  due  to  strong  X-electron  attenuation  in  the  GaAs  cap  layer.7 
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Fig.  11.  SD-BEEM  spectra  for  two  GaAs/GaInP2  (ordered  and  disordered)  samples  with  50  A  GaAs  cap  layers,  taken  at 
T  =  77  K.  For  clarity,  the  spectra  are  shifted  along  the  vertical  axis.  The  MSIS  model  calculations  are  also  presented 
(dashed  line).  The  top  inset  shows  the  photoluminescence  spectra  for  both  ordered  and  disordered  GalnP  samples, 
whereas  two  bottom  insets  show  the  schematic  of  the  GalnP  band  structure  (from  M.  Kozhevnikov,  V.  Narayanamurti, 
A.  Mascarenhas,  Y.  Zhang,  J.M.  Olson,  and  D.  L.  Smith,  Appl.  Phys.  Lett  75,  1 128  (1999)). 

The  SD-BEEM  spectrum  of  the  ordered  GaInP2  sample  presented  in  Fig.  1 1  show  a  very  important 
difference  from  that  of  the  disordered  GaInP2.  Namely,  we  observe  two  high-energy  peaks  instead 
of  one  peak  in  disordered  GaInP2.  We  assign  both  high-energy  peaks  to  be  associated  with  the  L 
valley  contribution.  In  CuPt-type  ordered  GaInP2,  one  of  the  four  L  valleys  folds  onto  the  T  point 
(hereafter  Lr),  and  the  other  three  are  folded  onto  the  D  point  (hereafter  Ld).  A  strong  repulsion 
between  the  r  valley  and  the  folded  Lr  valley  results  in  the  bandgap  reduction  and  in  the  increase 
of  the  F-  Lr  separation,  while  the  energetic  position  of  the  Ld  remains  almost  the  same.  Therefore, 
we  conclude  that  the  two  observed  high-energy  peaks  for  the  ordered  sample  are  the  contribution  of 
the  L  valleys  that  are  split  due  to  ordering. 

The  “folded”  zone-edge  bands  were  observed  experimentally  in  electro-reflectance  and  Raman 
spectroscopy24  measurements.  In  the  recently  reported  electro-absorption  experiments  on  ordered 
GaInP2  (r|  =  0.45)25,  an  additional  feature  was  observed  at  ~  0.48  eV  above  the  fundamental 
bandgap  transition,  and  this  feature  was  attributed  to  the  back  folded  L  conduction  band.  As  the 
ordering  decreases,  this  peak  shifts  to  lower  energies,  with  an  asymptotic  value  of  ~  0.33  eV  above 
the  fundamental  bandgap  transition  for  a  totally  disordered  sample.  In  Refs.23'25,  due  to  the 
selection  rules,  only  the  contribution  from  the  L  valley  folded  onto  the  T  point  was  observed.  In 
contrast,  we  observe  the  contribution  from  all  L  valleys  and,  as  a  consequence,  can  measure  directly 
the  F-L  separation  in  disordered  GaInP2  as  well  as  the  ordering-induced  L  valley  splitting  in  ordered 
GaInP2.  According  to  our  results,  A(F-L)  =  0.35  eV  for  a  disordered  sample,  A(F-Lr)  =  0.47  eV  and 
J(Lr-Lo)  =  0.13  eV  for  an  ordered  sample.  These  results  are  in  a  good  agreement  with  the 
theoretical  predictions.  Indeed,  as  pointed  out  by  Zunger, 26  it  is  possible  to  obtain  the  dependencies 
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A(r~Lr)  =  A(r-L\=o  +  0.5r|2  and  A(F-L)  n=o=0.35  from  our  experiment,  A(F-Lr)  =  0.475  eV  and 
A(Lr-Lo)  =  0.11  eV. 

This  then,  is  an  example  of  obtaining  entirely  new  information  on  band  structure  of  a  new  material, 
which  heretofore  was  not  possible. 


V.  GaSb  Quantum  Dots  Grown  on  GaAs 


GaSb  self  assembled  quantum  dots  grown  by  MBE  on  GaAs  exhibit  a  staggered  (type  II)  band 
lineup  with  a  potential  barrier  in  the  conduction  band.28'30  Therefore,  electron-hole  recombination 
in  these  structures  is  spatially  indirect,  so  that  optical  measurements  do  not  provide  adequate 
information  about  the  conduction  band  offset  between  the  GaSb  dots  and  GaAs.  Along  with  that, 
traditional  transport  methods  cannot  measure  this  local  band  offset  because  of  the  small  (—50  nm) 
lateral  dot  size,  and  any  offset  measurement  would  average  over  areas  with  and  without  dots.  In 
BEEM,  however,  carriers  are  locally  injected  into  semiconductor  structures  to  spectroscopically 
probe  buried  interfaces  on  a  nanometer  scale. 

Figure  12  shows  STM  and  BEEM  images  of  a  single,  GaAs  capped,  GaSb  dot.  In  the  STM  image, 
a  roughly  circular  feature  ~  50  nm  in  diameter  and  ~  5  nm  tall,  marks  the  lateral  position  of  the 
buried  dot.  The  area  in  the  BEEM  image  aligned  with  the  dot  profile  in  STM  is  darker  than  the 
surrounding  region,  implying  that  the  BEEM  current  through  the  dot  is  reduced  due  to  electrons 
reflection  off  the  dot’s  potential  barrier.  The  height  of  this  barrier,  i.e.  the  local  band  offset,  can  be 
extracted  from  the  changes  in  BEEM  spectra  between  the  on  and  off  cases.  The  on  dot  and  off  dot 
BEEM  spectra  of  several  dots  were  fitted  by  using  a  modified  Bell-Kaiser  planar  tunneling  model,31 
giving  a  local  conduction  band  offset  for  GaSb  dots  on  GaAs  of  0.08±0.02  eV.  For  details  of  this 
study,  see  Ref.32. 


Band  profile 


Au  film 


Images  BEEM  spectrum 


Tip  Bias  (V) 


Fig.  12.  (from  left  to  right)  Schematic  of  band  profile,  STM  and  BEEM  inages  and  characteristic  BEEM  spectra  for 
GaSb/GaAs  self-assembled  single  quantum  dots  (from  M.E.  Rubin,  H.-R.  Blank,  M.A.  Chin,  H.  Kroemer,  V. 
Narayanamurti,  Physica  E2,  682  (1998)). 


16 


VI.  GaN 


Although  GaN  exhibits  high-efficiency  optical  performance  despite  its  extremely  high  dislocation 
densities,33  the  dislocations  are  very  likely  to  have  significant  electrical  and  structural  impact  on 
high-power  electronic  devices. 


Fig.  13.  Collector  current  image  (left)  and  simultaneously  taken  STM  image  (right)  of  the  70A  Au  cap  layer  of  a  2500A 
x  2500A  area  (IT  =  InA,  VT  =  -1.8V).  The  small  high-current  area  (~500pA, bright  spot)  surrounded  by  areas  with  zero 
collector  current  (<  0.2pA)  shown  in  (a)  is  typical  for  the  investigated  GaN-films.  The  Au  layer  (right)  exhibits  a  pit 
(dark)  at  the  location  of  the  high-current  area  (bright  spot  (left))  (from  E.  Brazel,  M.A.  Chin,  and  V.  Narayanamurti, 
Appl.  Phys.  Lett.  74  (16),  2367  (1999)). 


The  recent  study  of  Brazel  et  al 34  has  revealed  high  BEEM  current  densities  and  low  effective 
Schottky  barrier  heights  in  small  localized  areas  around  dislocations  with  a  screw  component.  The 
GaN-films  used  in  that  study  were  grown  by  metal  organic  vapor  deposition  (MOCVD)  at 
atmospheric  pressure.  For  these  measurements  unintentionally  doped,  n-type  GaN-films  of 
thicknesses  ranging  from  1/xm  up  to  3/xm  were  deposited  on  c-plane  sapphire  substrates.  Fig.  13 
(left)  shows  a  typical  collector  current  distribution  for  the  MOCVD-grown  GaN  films  deposited  on 
c-plane  sapphire  substrates.  From  most  parts  of  the  sample  no  collector  current  was  detectable 
(dark  area),  and  only  in  small  localized  spots,  like  the  shown  light  area,  high  collector  currents 
could  be  observed.  In  this  representative  case,  the  measured  current  was  ~500  pA,  which  is  half  of 
the  injected  InA  tunneling  current.  To  localize  the  high  current  areas  on  the  sample  surface  and 
possibly  correlate  them  with  surface  features,  images  of  the  current  distribution  and  the  topography 
were  taken  simultaneously.  Fig.  13  (right)  shows  the  STM  image  of  the  Au  cap-layer,  which 
corresponds  to  the  collector  current  image  in  Fig.  13  (left).  The  70  A  thick  Au-layer  exhibits  a  -500 
A  wide  pit  at  exactly  the  same  position  of  the  high  current  area  in  the  collector  current  image.  An 
AFM-image  taken  of  the  same  sample  shows  several  of  these  pits  that  represents  the  surface 
terminations  of  mixed  and  pure  screw  dislocations.  Thus,  it  was  concluded  that  the  high  BEEM 
current  areas  are  located  at  dislocations  with  screw  character.  To  get  information  about  the 
Schottky  barrier  heights,  the  BEEM  spectroscopy  in  both  high  and  low  currents  areas  was  carried 
out.  The  low  current  areas  do  not  show  any  measurable  collector  current,  whereas  for  the  high 
current  areas,  we  have  found  barrier  heights  in  the  range  from  -0.3  eV  up  to  -0.95  eV. 
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Fig.  14.  Ic -Vt  spectra  for  both,  electron  and  hole  injection,  respecively.  Log-plots  like  the  one  shown  in  the  inset  were 
used  to  approximately  determine  the  threshold  voltage  necessary  for  current  injection  into  the  high-current  areas  (from 
E.  Brazel,  M.A.  Chin,  and  V.  Narayanamurti,  Appl.  Phys.  Lett.  74  (16),  2367  (1999)). 


The  threshold  values  vary  within  one  high-current  area  and  from  area  to  area.  In  addition,  it  was 
possible  to  reverse  the  tunneling  voltage  and  inject  holes  into  the  very  same  high  current  areas  as 
well,  with  the  observed  threshold  energy  of  ~0.3  eV  to  1.3  eV,  as  shown  in  Fig.  14.  The  observed 
thresholds  for  carrier  injection  do  not  seem  to  be  typical  Schottky  barriers,  but  rather  capturing 
barriers  of  trap  states  in  the  vicinity  of  the  dislocations.  These  results  definitely  suggest  that 
acceptor-  and  donor-like  trap  states  coexist  in  the  same  area,  and  are  being  filled  during  carrier 
injection  into  the  GaN.  The  maximum  threshold  energies  observed  for  electrons  (0.95  eV)  and 
holes  (1.3  eV)  add  up  to  2.25  eV,  remarkably  close  to  the  parasitic  yellow  luminescence  commonly 
observed  in  GaN.  In  the  future,  we  plan  to  conduct  more  experiments  to  answer  the  question 
whether  the  same  electronic  states,  which  are  causing  the  high  current  densities  are  also  responsible 
for  the  yellow  luminescence. 


VII.  Band  Offsets  in  Other  Systems 

Recently,  BEEM  spectroscopy  has  been  successfully  applied  to  determine  the  conduction  band 
offset  between  AlInGaAs  strained  layer  (active  layer  in  850  nm  vertical  cavity  lasers)  and  a  AlGaAs 
barrier  located  beneath  the  surface.35  Also,  BEEM  spectroscopy  has  been  used  to  study  both  the 
conduction  and  valence  band  offsets  between  InAs  and  AlAsSb  at  T  =  77  K  and  T  =  300  K.36  It 
was  found  that  with  the  addition  of  As  to  AlSb,  the  conduction  band  offset  between  InAs  and 
AlAsSb  decreases  despite  the  increase  in  the  bandgap.  The  resulting  decrease  in  the  valence  band 
position  causes  the  InAs/ AlAsSb  band  lineup  to  change  from  a  staggered  (type  II)  to  a  straddling 
(type  I). 
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SUMMARY 


Over  the  years  of  AFOSR  support,  we  have  successfully  developed  Ballistic  Electron  Emission 
Microscopy  (BEEM)  as  a  powerful  low  energy  electron  spectroscopy  technique  for  probing 
nanometer  size  electronic  objects  and  defects.  The  work  has  been  widely  published  and 
disseminated. 
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